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Abstract. From the hypothesis that at zero temperature the square root of the spectral continuum thresh- 
old so is linearly related to the QCD scale A we derive in the chiral limit and for temperatures considerably 
smaller than A scaling relations for the vacuum parts of the Gibbs averaged scalar operators contribut- 
ing to the thermal operator product expansion of the p° current-current correlator. The scaling with 
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A = y so CO /so (0), so being the T-dependent perturbative QCD continuum threshold in the spectral inte- 
gral, is simple for renormalization group invariant operators, and becomes nontrivial for a set of operators 
which mix and scale anomalously under a change of the renormalization point. In contrast to previous works 
on thermal QCD sum rules with this approach the gluon condensate exhibits a sizable T-dependence. The 
p -meson mass is found to rise slowly with temperature which coincides with the result found by means of 
a PCAC and current algebra analysis of the p° correlator. 
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1 Introduction 



QCD sum rules at finite temperature have been of intense 
interest since the pioneering work of Bochkarev and Sha- 
poshnikov [Q . Since thermal field theory lacks asymptoti- 
cally measurable states the empirical side, that is the spec- 
tral function of the QCD sum rule under consideration, 
must be parameterized. There are few guidelines for the 
parametrization of the spectral function by hadronic reso- 
nance and continuum contributions. However, as pointed 
out by several authors |§,|3[^] , the output of the sum rule 
at finite temperature or for nonvanishing chemical poten- 
tial depends on the hadronic model leading to a specific 
spectral function. 

Another complication due to finite temperature arises 
from the fact, that temperature is defined in a fixed ref- 
erence frame which hence is singled out. Therefore, the 
Poincare invariance of a field theory at T = is par- 
tially broken for finite temperatures. Residual 0(3) and 
translational invariance permit a wider set of operators 
to contribute to the operator product expansion (OPE) 
of the corresponding correlator. Using a background field 
method, Mallik g derived an OPE for the time-ordered 
thermal correlator of various quark field bilinears, where 
these new operators are included up to mass-dimension 
four. In contrast to the work of Hatsuda et. al. |2), which 
relies on the higher twist classification of operators origi- 
nating from the analysis of Deep Inelastic Scattering (DIS) 
and which for nonvanishing spatial momentum compo- 
nents allows for contributions of 0(3) non- invariant op- 



erators to the thermal OPE [J a systematic 0(3) invari- 
ant extension of the zero temperature OPE is obtained 
in Ref. . Thereby the authors consider the mixing of the 
new operators under a change of the renormalization scale. 
However, using the background field method, the Wilson- 
coefficients for the radiative corrections of mass-dimension 
six could not be calculated, and hence have been omitted 
in Ref. |J. These contributions are important, since they 
distinguish the vector from the axial vector channel and 
cancel a large part of the terms with mass-dimension four, 
yielding the experimental value of the p -meson mass at 
T = with good stability @ . 

In Refs. 0,§,§] the Gibbs average of the OPE is as- 
sumed to be saturated by the vacuum and the dilute ther- 
mal pion gas contributions. Thereby, both the vacuum and 
the pion states are taken to be temperature independent. 
This leads to a T-independent gluon condensate Q2j , which 
is in contrast to lattice measurements [|| and to results ob- 
tained in effective meson models with scalar glueball fields 
H|i"o|] . The main point of this paper is to explore the con- 
sequences of a T-dependent vacuum for the T evolution of 
the gluon condensate and the p -meson mass in the chiral 
limit. To this end, we parametrize the T-dependence of the 
vacuum implicitly through the effective spectral variable 
sq(T) which divides the hadronic part from the pertur- 
bative QCD domain of the spectral function. This ansatz 



1 At first sight this seems to violate rotational symmetry 
which is retained in the heat bath. Invariance with respect to 
0(3) transformations is, however, restored in Ref.|| when per- 
forming the phase space integrals over pionic matrix elements 
introduced through the definition of the thermal average. 
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can be justified for renormalization scales of the order of 
1 GeV and for temperatures considerably less than the 
fundamental QCD scale A (see Appendix A). 

The paper is organized as follows: In section II we list 
the basic results for thermal dispersion relations and for 
the spectral function in the p°-meson channel as already 
established in Refs. [|| and ||. Section III contains the 
thermal OPE for the invariant longitudinal amplitude of 
the p°-meson correlator, and we briefly discuss the behav- 
ior of the nonscalar contributions of mass-dimension four 
under renormalization. The treatment of thermal opera- 
tor averages is performed in section IV. With an implicit 
dependence on temperature, we derive a scaling relation 
with respect to A = y/ sq(T) / s (0) for T-dependent vac- 
uum averages of scalar operators. This relation is easily 
implemented for renormalization group invariants (RGI), 
and becomes quite involved if one has to regard a set of op- 
erators which mix and scale anomalously under a change 
of the renormalization point. In section V we write out the 
Borel transformed sum rule and, by performing a logarith- 
mic derivative with respect to the inverse squared Borel 
mass, obtain a sum rule for the p -meson mass which is 
the ratio of two moments. A numerical evaluation of the 
thermal sum rule is performed in section VI. Section VII 
summarizes and compares the results with those of previ- 
ous approaches. 



2 Thermal Dispersion relations, kinematic 
Invariants, and Spectral Function 

We consider the thermal correlator of the time-ordered 
(T) product of two vector currents in the /3°-channel 

T^{q,T) = Z~H I d^x e^Tr e~ l3H Tj ll (x)j v (0) , (1) 

where 



3n( x ) = o ( u 7m u - dj^d) and Z = Tr 



(2) 



Here H denotes the QCD Hamiltonian, and (3 stands for 
the inverse of the temperature T. A sum rule for this cor- 
relator can be derived as [pi 



V(Q ,|q|,T) 

Ql 



i 



, 2 ImT^ q , q , 
dq ° 72 7^2 — t anh(/3 9o /2) 



where {g^ v ) = diag(l, — 1, — 1, — 1), and the tensors P^, 
are given by 



p _ _ <r - - 

q z q 

n _ St~ ~ r = SdL 



(5) 



By evaluating TTi = u^u^T^ and II2 = T£ in the rest 
frame of the heat bath (u^ = (1, 0, 0^0)), one can solve for 
the invariant amplitudes 7] and T t [|6| as 



1 1 / 

T = ^n 2 , T t = -- [n 1 + %n 2 

q 2 2 V T 



and in the limit q — > one obtains 

T t (g ,q = 0) = q 2 T,(gb,q = a) 



(G) 



(7) 



Since the sum rule of Eq. (g) holds for each component 
of (T^), it also holds for TTi and 77 2 . With Eq. @ one 
obtains sum rules for the invariant amplitudes Ti and Tt. 
For example, 



Ti(ql,q,T) 



*r* N *f>«$ , with 



q'° 2 + Q° 2 



Ni(q'°,q,T) = Tr^Im^ tanh(/3 g '°/2) 



(8) 



As already indicated in the introduction, thermal field 
theory is handicapped by the lack of asymptotically mea- 
surable states, and therefore the spectral function (de- 
fined as the numerator of the integrand of the right-hand 
side of Eq. (||), that is Ni(q'° ,q,T)), must be modelled. 
In the following so(T) denotes a temperature dependent 
threshold that divides the hadronic part of the spectrum 
from the perturbatively accessible QCD domain. It is sug- 
gested that in the hadronic region of integration, 

< q' a < sq, the correlator is saturated by the p°- 
resonance and the two-pion continuum. Thereby, the pions 
are assumed to be noninteracting. The hadronic contribu- 
tions jjf and to the current of Eq. (|J) are obtained 
by using the ficld-current-identity 



31 



m p f p pl, where f p (T = 0) = 153.5 MeV , (9) 



and by appealing to the SU(2) flavor symmetry structure 
resulting in 

(10) 



3l = e 3bc 7r%7r 



-<& 



(3) 



Thereby, the pionic current is obtained by constructing 
the Noether current of the SU(2)y symmetry of a low- 
energy chiral Goldstone-field theory which to lowest order 
in the derivatives contains noninteracting fields. Since we 
only consider a dilute pion gas (one-pion states in Gibbs 
average) this free pion theory is relevant for our purposes. 

Employing the noninteracting, finite temperature (real- 
time) p°-meson and pion propagators |ll||l^,|l3| for the 
calculation of the unitarity cuts for the tree diagram of 
T^(q,T) = Q^Ti(q 2 ,oj,T) + Pp, u T t (q 2 ,uj,T) , (4) the p° contribution and for the thermal pion loop, the 



The frame of reference where temperature is defined moves 
with four- velocity u^. With this additional covariant one 
can, for example, define the Lorentz scalars uj = u^q 11 

and q = \/u> 2 — q 2 . Imposing current conservation and 
symmetry under exchange of p and v, the correlator of 
Eq. (pi) can be decomposed as pllrjl] 
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hadronic part of the spectral function has been calculated 
in Ref. ||. In contrast to Ref.[|| and following Ref.B we 
approximate the imaginary part of the correlator above 
the threshold s (T) by perturbative QCD (pQCD). For 
this perturbative piece thermal contributions are omitted 
since the fermionic distribution function rip remains small 
for the relevant temperature and energy range Q. In the 
limit q — > and after a Borel transformation in Qq with 
Borel mass M [Q the spectral side of the sum rule of Eq. 
(||) reads 

T t (M, T) = —2 / ds Im7](s, T) e" s '< IM " W(Vi / \2T)) 
withs = g '° 2 , (11) 



expansion for the invariant 7/ 



and 



48tt 2 J im 2 8tt 2 V 7T 

J t= / dse ~ 3/M ' 2 ^ / ds 
Js (T) 24tH 7 4m 2 

(e-'/«V + «(3 - « 2 )/2) n B (V5/2, T) , (12; 

where p — 1 GcV. Here the function v is dehncd as 
u(s,TO7r) = y/1 — 4m^/s, Ub denotes the Bose distribu- 
tion, and Jqfrpy Jq 9 are the vacuum (thermal) parts of 
the spectral integrals due to the irir and pQCD continua, 
respectively. 



3 Thermal Operator Product Expansion 

In this work we use a thermal OPE for the invariant ampli- 
tude T\ which combines the results of Ref. [|| and Ref. Q . 
In Ref. H the expansion is only carried out up to operators 
of mass-dimension four. The nonscalar 0(3) invariant con- 
tributions are expressed in terms of diagonal combinations 
with respect to the anomalous mixing matrix, resulting 
in a renormalization group invariant (RGI) (that is the 
total energy density) and a renormalization group non- 
invariant (RGNI) contribution. The drawback of the OPE 
of Ref.|| lies in the fact, that mass-dimension six con- 
tributions are omitted. These terms are important since 
they cancel a large part of the nonperturbative correction 
of mass-dimension four yielding at T — the experimen- 
tally measured p -meson mass. In addition, it is the mass- 
dimension six part of the OPE that distinguishes at T = 
the vector from the axial vector channel |7]]. Therefore, we 
use together with the standard T = scalar operators the 
results of Ref. Q for the nonscalar mass-dimension six con- 
tribution. In the chiral limit, for q — > 0, and for two light 
quark flavors (u and d) 0, we then obtain the following 

2 In order to compare our results with the lattice data of Ref. 
H for two quark flavors we constrain ourselves to QCD with 
two light quark flavors throughout the paper. 



T l ( q ° 2 ,c l = o,T) 
1 
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87r* \ p 
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IT 



7T 



1 ,0.,, 



F%F^ a ) T {p) + 



7T 
7T 

8ni 
3Q 



6 x a s {wy^t a u - d llll5 t a d) 2 ) T {Q ) 



-^{a^u-y^u + d lt ,t a d) 2 ) T {Qo) + 



6 {wjoDqDqDou + d-y D DoDod}T{lJ,), 



(13) 



where the nonscalar operators are understood to be sym- 
metrized and to be made traceless with respect to the 
Lorentz indices. The generators t a of color SU(3) in the 
fundamental representation are normalized to Trt a t b = 
2 6 ab . There are, in principle, also 0(3) mixed quark-gluon 
operators of twist 4. However, their there is no experi- 
mental clue about their matrix elements as pointed out in 
Ref. g. There it was also indicated that bag model esti- 
mates of the nucleon matrix elements of these operators 
yield very small values and that the treatment of pions in 
the bag model is rather questionable due to their collective 
nature. Following Ref. |2j we simply omit these operators 
to obtain Eq. @. 

The Gibbs average in Eq. (|l^) is approximated by the 
vacuum and the dilute pion gas contributions in Ref. [Q. 
As the authors point out the scalar part of the Gibbs 
averaged OPE of Eq. ( fi"3| ) then has the same structure as 
the result for the vector correlator obtained in Ref. EEJ. 
Based on PCAC, current algebra, and the LSZ reduction 
formula the authors of Ref. |15 find to order T 2 a mixing 
of the correlator of the vector with that of the axial vector 
channel due to finite temperature. 

To one loop and for two quark flavors the constants b 
and S are given by [Em 



6= A- 

3 V 3 



2), 



11 



(14) 



At fj, = 1 GeV we take a s (p 2 ) = 0.36 f|. Note, that 
in the case of mass-dimension four the operators a priori 
renormalized at Q° are already expressed by operators 
evaluated at p = 1 GeV. Only for (-y#oo — ®oo)t does 
this process of rescaling generate a renormalization group 
logarithm due to the nonvanishing anomalous dimension 
5 of this operator. 

The operator ^F£ V F» V ' a is an RGI. In the chiral limit, 
the fermionic and gluonic parts of the 00-component of 
the traceless energy- momentum tensor {0^ u ) of one quark 
flavor QCD §,fl7| are 

#oo = WJoDoq , 



^oo 



'^OA^Oa 



^ffOO F£ X F K a , 



D = d -igA^ 



(15) 
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with (9qo given by 



9oo - # 00 



(16) 



Since 9^ v is a conserved quantity the thermal average 
of #oo is also an RGI. The expression for the covariant 
derivative Do in Eq. ( |l5|) contains the zeroth component 
of the gauge potential A$ and the color SU(3) genera- 
tors t a in the fundamental representation normalized to 
Trt a tb = 2(L f) 

In Eq. (|13|) both mass-dimension six scalar operators 
are RGNI, and we will consider their anomalous scaling 
and mixing in the next section. As for the 0(3), non-scalar 
mass-dimension six contribution to Eq. ( |l3| ) we neglect the 
anomalous scaling of the corresponding operator g . It is 
known to belong to a set of operators which mix under 
renormalization. This set also contains mixed quark-gluon 
contributions, and following Ref. @J we will omit them. 



4 Gibbs averages 
4.1 Scalar operators 

As was suggested by several authors Jl[|,|| the Gibbs av- 
erages of the scalar operators of Eq. (R3 ) can be saturated 
by the vacuum and the one pion contributions. Thereby, 
both the vacuum and the pion states were assumed to ex- 
hibit no temperature dependence 0]. In the chiral limit 
the only quantities entering the sum rule, which can po- 
tentially describe thermal pion properties, are the pion de- 
cay constant and the pion matrix elements of twist two 
operators. As for the former the T dependence has been 
obtained in the imaginary time formulation of thermal 
chiral perturbation theory in Ref. jl8|. There, the result 
for fir(T) is a notably decreasing function of temperature 
for T > fx. On the contrary, the lattice simulation of 
Ref.|l9|] obtains a nearly T independent pion decay con- 
stant up to the critical temperature which is understood 
as an artefact due to the large pion mass (~ 400 MeV) 
used. As for the T-dependence of the twist two pion aver- 
ages no information is available so far, and we have to use 
the pion-in-vacuum parton distributions to estimate the 
corresponding matrix elements. In accord with previous 
sum rule investigations at finite temperature |^,g,|l| and 
for consistency we then have no choice but to assume a 
T-independent pion decay constant. 

The motivation for a T-dependent vacuum part in the 
Gibbs average of scalar operators emerges from a com- 
parison of lattice data for the gluon condensate || with 
the result of Ref. . The former approach indicates a siz- 
able decrease above a critical temperature of about T c — 
140 MeV, whereas in the latter case the gluon conden- 
sate exhibits practically no temperature dependence even 
when leaving the chiral limit (less than 0.5% decrease at 
T = 200 MeV)0. On a more phcnomcnological level, the 

3 On the lattice and at T = the ground state average 
of the local two gluon operator F^ V (Q)F£ 1 ' (0) in contrast to 
the gluon condensate used in an OPE contains also perturba- 



T-dependence of the gluon condensate has been obtained 
from the effective potential of theories based on the non- 
linear cr-model, where the mesons are coupled to a scalar 
glueball field to mimic the breaking of scale invariance 
by the QCD vacuum Depending on the normaliza- 

tion conditions used for the Bag constant and the glueball 
mass, in these calculations the T-evolution of the gluon 
condensate exhibits a strong decrease above critical tem- 
peratures ranging from T c = 140 — 400 MeV. 

The scalar operators appearing in the OPE of the ther- 
mal current correlator read 



dim 4 
dim 6 



TT 



Or, 



ira s (u^fjf'u + d^fy.f'd)' 



(17) 



Hereby, 0% is an RGI, and 0§ii CI, 2 can b e expanded 
into a basis of scalar four-quark operators Pi, P%, Pq 
[Ol with 



Pi 

P 2 
P 3 
Pi 
P> 
Pe 



4>Ll^L$Lln^L + {L -> R) 

TpLl^^L^Rl^^pR , 
(«L(fi),4(R)) T , 



(18) 



where t a and X b are the color SU(3) Gell-Mann and flavor 
SU(2) Pauli matrices, respectively. They are normalized to 
Trt a t b = Tr\ a X b = 2S ab . The left-(right-) handed spinors 
are given by 

^ L{R) = \{l± l5 )i) . (19) 



Using the relation 



r ij T mn ~ 2(Sin8j m 



8ij S mn ) 



(20) 



for the SU(N) generator matrices r c (Trr a r b = 28 ab ) and 
applying Fierz transformations, one obtains the following 
decomposition of the flavor singlet parts^] (for a derivation 
see the Appendix B) 



O s 61 = ira s - 



5 32 „ 
-Pi + — P 3 
3 9 



2 P 



O, 



6,2 



na s (P 4 + 2P 2 ) 



(21) 



tive contributions which must be subtracted for a direct com- 
parison. However, at finite temperature the perturbative part, 
which contributes to the partition function for momenta larger 
than, say, 1.5 GeV, is for T < 200 MeV severely Boltzmann 
suppressed and hence hardly influences the T dependence of 
the thermal average. 
4 Notationally not quite correct we also refer to them as 
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The sets of operators P\ , . . . , P4 and P5 , Pg mix indepen- 
dently under renormalization with the respective one-loop 
mixing matrices S and 5 |2C|j 



/ 3/2 
16/3 17/3 
-2/3 
\ -20/9 




7 16/3\ 
3/2 J 



. (22) 

The eigenvalues (proportional to the anomalous dimen- 
sions of the diagonal combinations) read 



£1 = 7.043 , 6 2 = 3.501 , S 3 
= 8 , and Sq = — 1 . 



-2.891 ,S 4 



4.097 



(23) 



In Appendix A we argue that for temperatures consid- 
erably smaller than the fundamental QCD scale A (A 
200 MeV HUl), in the chiral limit, and for a renor- 
malization scale Qo of the order of 1 GeV it should be 
possible to describe the T dependence of the thermal av- 
erage of a scalar operator implicitely through that of the 
T-dependent spectral continuum threshold sq(T). Meeting 
the above premises, we can derive a scaling relation with 
respect to A = \J 'so(T) / s (0) for the vacuum average of a 
given scalar and diagonal operator O with mass-dimension 
d and anomalous dimension one-loop coefficient So- In the 
spirit of a random phase approximation, where correla- 
tions between particle-hole excitations in the ground state 
are included, we set up the QCD vacuum state as an ex- 
pansion in terms of fc-particle (off-shell) quark-antiquark 
and gluonic fluctuations with vacuum quantum numbers, 
that is 



\0)=J2\k), with 

k 

\k)= Jd 4 Pijd 4 p 2 ---Jd< 



Pk 



,(Pi,Pa, -,Pk) \PV,h) ■ ■ ■ \Pk\ik) 



(24) 



Hereby, ij is a collective index labelling the particle species 
(fermion or boson), spatial quantum numbers, flavor (if 
fermionic), and color. The Ci lt i at .. i h denote the corre- 
sponding expansion coeffcients. Then the vacuum average 
of O has a representation of the following form 



(0|O|0)(Q ) = E E / / --Jd* 



Pk 



k »i,ia, •••>** 
fii,i 2 ,...,i k (Pl,P2, —,Pk',Qo) , 

(25) 

where the functions fi u i 2 ,...,i h have mass-dimension d— 4k, 
and Qo denotes the scale at which the operator O is renor- 
malized. By asymptotic freedom, the integrand of Eq. ( p5[ ) 
will be strongly suppressed, if the momenta p±, ...,pfe are 
hard. The finite value of (0|0|0)(Qo) has its origin in the 
strong dynamics of soft fluctuations in the vacuum [M . A 
criterion distinguishing soft from hard fluctuations should 



roughly be given by the scale sq. For sizable contributions 
to the integral of Eq. J25| ) we assume that the time- and 
space-like virtuality p 2 of each of the momenta p\, ...,pk 
be less than so and, in addition, that (p^) 2 be less then sq 
(recall that through the presence of the heat bath (p®) 2 
is formally raised to a Lorentz scalar). With the above 
premises and the results of Appendix A we may assume 
that temperature effectively acts on the vacuum only im- 
plicitely through so and that there is a linear relation be- 
tween the T dependent QCD scale At and \J sq(T). Thus 
there are only two independent scales to be considered: 
so(T) and Qo- Performing the integrations over the spa- 
tial components of the momenta in Eq. ( |25j ) we obtain 



(0|O|0)(Q 

e e rZ:<&~ 

k ii,i2,...,i* - V*°(0) 



\Ao(o) 



dpi 



foil ,22 ,. 


..,iM,P°2, 




({I' 




\p\ 


Pl \ 


\Qo' 


\/«o(o) J 



P? Pi 



Pl J 1 Qo v^oW 



Qo 



(26) 



In Eq. ( |26|) the functions /ii 1 ,i 2 ,...,i s . are homogeneous func- 
tions ofpi, (i = 1, • • • , k) with mass-dimension d— k, and 
9iiM,---,ik are dimensionless functions of their dimension- 
less arguments. The T-dependent vacuum average (not to 
be confused with the Gibbs average) is then given as 

(Q|O|0)r(Q ) = 



EE/ dpi- dpi 

hh,i 2 ,...,i k (Px,P2i — iPfc) x 



Pl Pl 



' rfi ri° 1 J rf> 

Pl Pk-1 I I Pl Pl 



.Pfe Pfc 



Qo ^{T)\ VW), 



With 



A = 



l so(T) 
so(0) ' 



(27) 



(28) 



and from comparison of Eq. ( Pq ) and Eq. (|27j ) one easily 
obtains 



(0|O|0) T (AQ ) = A d (0|O|0)(Q ) 



(29) 



Eq. (|2^) relates the average of a diagonal operator O with 
mass dimension d taken in a temperature destorted vac- 
uum and renormalized at AQo to the T = vacuum aver- 
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age renormalized at Qq. This scaling relation should em- 
body a good approximation for the thermal vacuum aver- 
age, provided that the above premises of massless quarks, 
T considerably smaller than A, and Qq of the order of 1 
GeV are fulfilled. 

For an RGI operator O (i5o=0) the rescaling from AQo 
to Qq is trivial (for example Of). The situation becomes 
more involved for the scalar operators Og l and C| 2 of 
mass-dimension six given in Eq. (|l7j), and we will focus 
on them now. 

The perturbative one-loop renormalization group rescal- 
ing from AQo to Qo for the diagonal operator O is given 
by 



_>/.- 



(O|0|O) r (Qo) = 



log(AQ ) 2 M 2 



2\ -(So/b) 



(O|0|O)t(AQ o ) 



log Ql/A 2 

_ (30) 
where b is defined in Eq. (|l4|). In the OPE, operators renor- 
malized at Qq are expressed by operators renormalized at 
a common reference point fi. Eq. (^9| ) implies that first we 
rescale from AQo to Qq which is accomplished by 



\O\0) T (Q a ) = X 



d /log(AQ ) 2 M 2 



-(Sa/b) 



\ logQ^M 2 
and then express (0|O|0)(Q ) by (0|O|0)(^) as 

'log(AQ ) 2 M 2X ~ (5o/b) 



<0|O|0)(Q ), 
(31) 



(0|0|0) T (Oo) 



log Q 2 /A 2 



log Ql/A^ So/b 



log n 2 /A 2 



\O\0)(ji) . (32) 



In order to apply the following identity for the Borel trans- 
formation Lm Q 



1 



1 

Ql 

(JL 

r{k) \m 2 

l 



i 



HQl/A 2 ) 



l 



ln(M 2 /yl 2 ) 



O 



1 



ln(M 2 /yl 2 ) 



{co + Cl (A 2 -l ) + C2 (A 2 -l) 2 + 

0((A 2 -l) 3 )}(ln( M 2 M 2 )) (1 - fe/b) (0|O|0)(m). 

With M 2 fa 0.75 GeV 2 and yl 2 = 0.04 GeV 2 the coeffi- 
cient C2 in the expansion of Eq. (34) is roughly given by 
l/2c\. So even for a value of A 2 as ow as 1/2 we obtain 
a suppression for the quadratic as compared to the linear 
term in Eq. ( |3~i| ) by a factor of four. 

In the case of the two sets Pi, P4 and Ps,Pe of Eq. 
(Eq) the mixing under renormalization results in 




\0)t(M) = TDT 1 (0|7ra s 
\0) T (M) = f DT-^OIttcks 




|o)W 
|o)(m) , 

(35) 



where the matrices T and T, transforming the matrices 5 
and 8 of Eq. (E2h , respectively, are given by 



0.1965 -0.0556 0.0496 0.6838 
0.9224 -0.1298 -0.0955 -0.5003 
-0.0008 0.4786 -0.5481 0.3554 
-0.3324 -0.8666 -0.8295 0.3947 



/ 8/9 -16/27\ 

W 6 8 / 9 / 



(36) 



The diagonal matrices D = diag(£>i , D2, D3, D4) and D 
diag(Z?5,I?6) have matrix elements of the form 

A = A 6 j(ln(M 2 M 2 )) (WH - 
(Si/b) (HM 2 /A 2 )f i/b - 2) (A 2 -l) + 



^(^ + 6(l + lnO, 2 M 2 ))) (ln(M 2 M 2 )) 

(A 2 -1) 2 + 0((A 2 - 
(i-Si/b) 



(Si/b-3) 



I) 3 )} 



(33) (hVM 2 )) 



(* = 1,.» ,6) , 



(37) 



sxpand (^) 



|C|0)t(Qo) about A 2 
order and perform afterwards the Borel transformation in- 
dicated in Eq. (p3|). The result is 



where 5i, ...,Sq are given in Eq. (|23]). Note that the one- 
1 up to quadratic^ 00 ? sca ling of a s is also included in D and D. When 



\2/,- 



i) U|tl|(i::,((J„l 



-LY 



r(k) \m 2 J 

(So/b) (ln(M 2 /A 2 )) 



{(MM 2 /A 2 )) 



So/b-l 



(A 2 - 1)+ 



^L[5 +b{l + H^/A 2 )}] (ln(M 2 M 2 )) 
0((A 2 



(So/b-3) 



l) 3 )}(ln( / , 2 M 2 )) (1 - fe/b) <0|O|0)( M ) 



evaluating the sum rule numerically we will consider a 
truncation of Di after terms of zeroth and second order in 
(A 2 — 1) to test the sensitivity of the T-evolution of m p , 
So, and the gluon condensate. 

Making use of the vacuum saturation hypothesis of 
Shifman, Vainshtein, and Zakharov (SVZ) [|14| 

(o|^ri#r 2 v|o) = N- 2 [Trr 1 Trr 2 - Tr/Y^KO^O) 2 , 

(38) 

where N = 4 x N f x N c , and the r t , (i = 1,2), are 
direct products of Dirac, flavor, and color matrices, one 
_caa_easily verify Q that 

(0|P 2 |0)( M ) = ^(0|Pi|0)(/x) , 
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(0|P 3 |0)( M ) = (0|P 4 |0)(/x) = 
(0|P 5 |0)(/i) = ^<0|P 6 |0)( M ) 



, 



(39) 



With Eqs. (H), (H), @, and © we obtain a descrip- 
tion for the thermal vacuum contributions of scalar oper- 
ators to the Gibbs averaged OPE of Eq. ([l3|). 

As already mentioned in the beginning of this subsec- 
tion and following Ref. 0] we restrict the thermal trace 
of Eq. (0) to the contribution of one-particle pion states 
(the spectral integral over the pion energy can savely be 
extended to infinity because of the strong Boltzmann sup- 
pression of the integrand). We use the results of Ref. [g| 
which for 0\ were obtained by appealing to the trace 
anomaly of the QCD energy-momentum-tensor. In the chi- 
ral limit there is no explicit temperature dependence of 
the thermal average of the two-gluon operator 0% due to 
pion contributions. Applying the soft pion theorem twice 
and using the vacuum saturation hypothesis of Eq. ( |38| ) , 
according to Hatsuda et al. Q one obtains for the pionic 
part of the Gibbs average of the operators 0% 1 and 0% 2 



E 

a=l 

128 
~8T 



<""<">'^-. 



na s {qq) 



1 



3 T< 



:06, 2 K(P)} n B (\p\/T) 



(40) 



In the chiral limit we have/,, = 88MeV§], n,B is the Bose 
distribution, and q indicates a single light flavor quark 
field. As in Ref.Q we only consider perfect vacuum satura- 
tion since for the finite temperature evaluation we are only 
interested in relative changes as compared to the T = 
case. 



4.2 0(3) invariant Operators 

The 0(3) invariants of the OPE of Eq. ( |l3| ) have vanishing 
vacuum expectation values and therefore we only have to 
consider their pionic matrix elements. In the chiral limit 
we encountered the following 0(3) invariant operators of 
mass-dimension four in the OPE of Eq. (O) 



dim4 : O 



0(3) 



7 m> 



The operators 



0(3) _ 16,/ 

,0(3) 



nil 
7 (H) 



(41) 



and 4 2 are diagonal combinations 

and 

,,n n*|. The pion matrix elements of the quark part Nf0Q O 



3/ 



with respect to the anomalous mixing matrix of 9 00 

R3 



(Nf denotes the number of quark flavors considered) and 
the gluon contribution 9q Q to the total energy density 9oo 
are estimated to be equal and can be calculated from the 
valence parton distribution in the pion as found by Glueck 



5 For diagonal operators, we can, by means of the scaling 
relation of Eq. (^9|) , express their T-dependent vacuum average 
by a scaling factor times the T = vacuum average which 
vanishes due the boost invariance of the vacuum state. 



et al. Hi- Thereby a low-energy (Q 2 = 0.25 GeV 2 ) va- 
lence like parton distribution fitted to experimental data 
(direct photon production and Drell-Yan) is evolved at 
one and two loop order to the sum rule scale of about 
Q 2 = 1.0 GeV 2 . One obtains 



E/ ^ (n a (p)\e f 00 K(p))(Q) n B (\p\/T) 

X>2 ( " +d) (Q), 



TT 2 T 



120 



(42) 



with 



.4 



a(u-\-d) 



(lGeV) = 0.972 , Va 



(43) 



At mass-dimension six we consider according to Ref. | 
only the following twist two 0(3) operator 



dim6 



0° {3) = i (u l0 D Q D D Q u +(u^ d)) . (44) 



Again with the model of Glueck et al. and in the chiral 
limit the result for the pionic contribution to the Gibbs 
average of O?^ has been determined in Ref. 0] as 



3 


f d 3 p 
' 2|p| 








. 63 



>))(Q) n B (\p\/T) 



J2 A 



a(u-\-d) 



(Q) 



where 



.4' 



a(u-\-d) 



(1 GeV) = 0.255 , Va 



The pion matrix elements of the operators O 



0(3) 
4,1 i 



o 



(45) 
(46) 

0(3) 
4,2 ) 



and 0(3) 6 are of order T or higher which led the authors 
of Ref. [2] to omit them in the sum rule since otherwise 
higher orders in T would also have to be calculated for the 
scalar contributions which a priori are of order T 2 . We do 
not agree with this power counting argument. The fun- 
damental approximation for the inclusion of pions in the 
Gibbs average is that of a dilute pion gas. Higher orders 
in T would come in by considering interactions between 
pions as it was shown in Refs. |l^,^]. Since the pionic 
matrix elements of the above nonscalar operators happen 
to be of order T 4 and higher already when interactions 
are switched off the ommission of these operators would 
have the effect of artificially introducing pionic interac- 
tions. Hence we do not omit the nonscalar operators in 
the sum rule. 

We now have all the ingredients to the Gibbs averaged 
OPE of Eq. @. 



5 Sum Rule 

Performing a Borel transformation of the OPE of Eq. ([l3|), 
we obtain in the chiral limit the thermal sum rule for the 
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7KM) = JL(i 



2 

TT 
1 



(Ooo)t(p) + 



a s (AP 



-5/6 



1 



5 32 
, -P 4 + — Ps - 2 P 5 ) T (M) + 
M 6 Vl2 3 9 ; v ' 

I 471"'. 

— (P 4 + 2 ft)r(M) + — (q-foD D Doq)T{v) 



(47) 



where the dispersion integral for TJ(Af) appearing on the 
left-hand side of Eq. @ is given by Eqs. (|l2|) and (|ll|). 
The scalar mass-dimension six operators Pi, P3, P4, Pg renor- 
mahzed at M can be expressed by the operators P2 and 
P5 renormahzed at /z = 1 GeV with the help of Eqs. ( |35|) 
and (§|). 

In order to eliminate the coupling f p we solve Eq. ( |47| ) 
for the p -meson term 

p ee m y m > 12 , 

substitute r = 1/Af 2 , and perform the logarithmic deriva- 
tive of R [E3lE41 to obtain the ratio of moments 



d 

mJT, s , t) = -^-log P(T, s , r) 



dr 



(48) 



which we will consider. 



6 Numerical evaluation 

In this section we discuss the numerical evaluation of the 
sum rule of Eq. (48). Thereby, we employ the following 
values for the one-flavor quark and the gluon condensate 



■lw|o) (P 



(p 



1 GeV) = -(250 MeV) 3 
1 GeV) = 0.012 GeV 4 . 



(49) 



As a result of vacuum sum rule calculations in a variety 
of channels the central value for the gluon condensate is 
actually believed to be twice as high as the value stated in 
Eq. ( |49| ) with an error of about 50 %. The central value for 
the quark condensate indicated here is rather up to date 
and exhibits a smaller error of about 30 % [M. However, 
the higher values for the gluon condensate have been ob- 
tained using numerical correction factors n of about k = 2 
which multiply the four-quark condensate derived from 
the vacuum saturation hypothesis |14|l. For exact vacuum 
saturation (k = 1) the values of Eq. (p9[) for the quark and 
gluon condensates nicely reproduce the p -meson mass in 
the vacuum (m p {T = 0) « 760 MeV). Since we are only 
interested in relative changes of the spectral parameters 
and condensates for T > as compared to the zero tem- 
perature case we will stick to the values of Eq. (|49|) j|] . 
Our strategy in evaluating the sum rule of Eq. (081) is as 



follows: 

Since r is not a direct physical observable, it is chosen as 
the stationary point t s of m 2 p {Q, So(0), t) for a given So at 
T = 0. 

For a value of So(0) = 1.5 GeV 2 an additional contri- 
bution to the spectral function due to a\ — n production 
is regarded as a part of the pQCD continuum. This is rea- 
sonable as long as the Borel parameter r is larger than 
1 GeV -2 since structures stemming from the a\ — it pro- 
duction and possible radial excitations of the p -meson are 
then sufficiently suppressed in the spectral integral [^6| . 

There is a pronounced minimum for a numerically ob- 
tained value of A/ 2 = l/r s ss 0.73 GeV 2 corresponding to 
a p -meson mass of 745 MeV in the chiral limit considered 
.here. 

For finite temperatures T we determine the value of 
so(T) from the stationary point of m 2 (T, so(T), t) at the 
same value r = r s as in the zero temperature case. Per- 
forming this calculation at a sufficient number of T points 
yields the temperature evolution of m 2 , so and, with the 
parametrization of Eq. ( |29| ) , also the temperature depen- 
dence of the gluon condensate. The usefulness of the pro- 
cedure to determine the T-evolution for a fixed r s was 
checked by applying this technique to the p° sum rule 
given in Eq. (4.4) of Ref. § (s (0) = 1.5 GeV 2 ). Our 
method yields the same results as found in Ref. ^ , where 
a more elaborate evaluation analysis (averaging over a In- 
dependent Borel window) was used. 

We consider the following two cases: (a) naive rescaling 
[] and (b) renormalization group rescaling when including 
terms up to quadratic order in A 2 — 1 in the expansion of 
Di given in Eq. (|37]). 

The temperature evolution is determined in steps of 
AT = 3 MeV. For s (0) = 1.5 GeV 2 and in the case {b) 
the p -meson mass exhibits an increase of about 17.5% 
from its zero temperature value at a 'critical temperature' 
of T c — 157 MeV which is in contrast to the result of Ref. 
||. Thereby, T c is obtained by demanding the following: 
If a variation of s$(T) with respect to so(T— AT) is larger 
than 0.3 GeV 2 to produce the minimum of the Borel curve 
for m p at r s , we set T = T c . This critical value merely indi- 
cates the breakdown of our sum rule analysis and may not 
coincide with the critical temperature of a deconfinemcnt 
phase transition. There is no strong dependence of T c on 
the initial value of so for the case (6) - T c (so(0) = 1.2 
GeV 2 )= 157 MeV, T c (s (0) = 1.5 GcV 2 )=157 MeV , and 
T c (s (0) = 1.8 GeV 2 )=163 MeV. Fig. 1 shows the re- 
sults for m p as a function of the temperature, where for 
s (0) = 1.2 GeV 2 and s (0) = 1.8 GeV 2 only the case (6) 
has been considered. There is practically no difference in 
the evolution of m p when including the quadratic terms 
in A 2 — 1 for the Di of Eq. (|37|) as compared to the linear 
and to the zeroth order approximation. The deviation at 
T = 160 MeV is then at most 5 MeV. 



6 Naive rescaling means, that, according to Eq. (^), the 
anomalous renormalization group rescaling of the operators 
(powers of logarithms) are suppressed, leaving only the factor 
A 6 at mass-dimension six. 
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Fig. 2 indicates the temperature dependence of the 
pQCD threshold so(T) for the three initial values and 
the case (b). In addition, we show the case (a) choosing 
so(0) = 1.5 GeV 2 . Again, there is practically no difference 
between the results for truncations of Di in zeroth, linear, 
and quadratic order. Hence the logarithmic corrections in 
the Borel parameter r seem to have a much greater ef- 
fect on the T-evolution of m p and sq than the logarithmic 
corrections in A 2 . For so(0) = 1.5 GeV 2 the behavior in 
the case (b) is qualitatively similar to that found in Ref. 

However, the influence of sq(T) on the sum rule is 
quite different in our work. Besides the truncation of the 
hadronic part of the spectral function, s (T) also scales 
the vacuum part of the Gibbs averages in the OPE. In the 
case (a) we obtain an increase of sq (T) and hence a rise of 
the gluon condensate with increasing temperature which 
is in contradiction to the results obtained on the lattice 
and by effective meson models |||^,|To|j. 

In Fig. 3 we finally show the temperature evolution of 
the gluon condensate normalized to its T = value which 
due to the scaling relation of Eq. (p9f), the invariance un- 
der a change of the renormalization point, and the explicit 
T- independence in the chiral limit (according to Eq. (p9|)) 
is proportional to so(T) 2 . We indicate the same combina- 
tions of initial values sq(0) and cases (a) and (b) as for 
the T-evolution of sq. Fig. 3 indicates a universality of the 
T evolution in a sense, that for different values of so(0) 
almost identical results are obtained. 



7 Summary and Discussion 

The main concern of this paper was the investigation of 
the consequences of a temperature dependent vacuum for 
the T-evolution of the p -meson mass and the gluon con- 
densate in the chiral limit. Thereby, we used the method of 
thermal QCD sum rules. The results obtained in Refs. [^|, 
|(| for the thermal Operator Product Expansion containing 
also 0(3) invariant contributions and for the thermal p° 
spectral function have been combined in our calculations. 
Following Ref. the Gibbs averages of local, gauge in- 
variant operators contributing to the thermal OPE of the 
p° current-current correlator were saturated by vacuum 
and one-pion matrix elements. For lack of better knowl- 
edge as far as the T-dependence of pionic twist two matrix 
elements is concerned wc had to work with the T — par- 
ton distributions of Ref. |^l| . Consistency then demanded 
the use of a T-independent f n , although the T-dependence 
of this quantity is known for low temperatures from chiral 
perturbation theory . 

In contrast to previous sum rule calculations we sug- 
gested to account for a temperature dependence of the 
vacuum part of the Gibbs average. This was motivated by 
the observation, that the gluon condensate remains prac- 
tically T-independent up to temperatures of 200 MeV in 
the sum rule analysis of Ref. [Q , while a lattice measure- 
ment yields a drastic decrease at a critical temperature T c 
of about T c « 140 MeV || . To resolve this contradiction, 
we derived a scaling relation for the thermal vacuum av- 
erage of scalar operators with the T-dependent spectral 



pQCD continuum threshold sq. Thereby, the chiral limit, 
temperatures considerably smaller than the QCD scale A, 
and a renormalization point of the order of 1 GeV were as- 
sumed. The implementation of the above scaling relation 
is trivial for renormalization group invariant operators and 
becomes rather involved for operators which mix and scale 
anomalously under a change of the renormalization point. 
On the spectral side we used a narrow width approxima- 
tion for the p° resonance since including a finite width 
r leads to ambiguities in the determination of m p and 
r as was shown in Ref. || for the case of finite density. 
One would expect, that this is also true for finite temper- 
atures. On the other hand, the effect of the scaling of the 
vacuum averages can only be isolated if one compares the 
results with previous calculations also using the narrow 
width approximation. In analogy to a previous sum rule 
calculation in the p° -channel (see Ref. M ) we employed a 
thermal pion continuum for the hadronic part, while the 
high energy tail was approximated by a T-independent 
pQCD continuum. 

As compared to the vacuum case our sum rule calcu- 
lations indicate a rise of the p -meson mass (17.5% at T c ) 
and a decrease of sq and of the gluon condensate with 
temperature. Thereby, the critical temperature T c , where 
the sum rule analysis breaks down, is T c w 160 MeV. This 
"critical" temperature is rather close to the QCD scale A, 
and therefore we do not expect the scaling relation of Eq. 
( p9|) to be a good approximation for the thermal vacuum 
average of scalar operators. The increase of the p -meson 
mass contradicts the sum rule results obtained in Refs. jl], 
||,^3|, where a monotonic decrease of this quantity is ob- 
tained. In the difference sum rule analysis of Ref. Q the 
p -meson mass was also found to increase slightly (Am p w 
5 MeV) up to a temperature of about 125 MeV while de- 
creasing thereafter. There are indications for a slight rise 
of the p -meson mass up to T = 150 MeV (Am p ss 30 
MeV at T=150 MeV) from a microscopical calculation of 
the spectral function in the p° channel using an effective 
p — it Lagrangian |27| . We emphasize again, that the crit- 
ical temperature of our calculation does only indicate the 
breakdown of the sum rule analysis based on a scaling rela- 
tion for the scalar condensates and may not coincide with 
the critical value for a deconfinement phase transition. It 
is interesting to note, that a calculation in the Nambu- 
Jona-Lasinio model yields at T ss 160 MeV and already 
at nuclear saturation density a relative increase of the dy- 
namically generated constituent light-quark mass which is 
of the same order as that of the p -meson mass in our cal- 
culation (H,^. In Ref. jl5|], where by means of PCAC, 
current algebra, and the LSZ reduction formula a mixing 
of the vector with the axial vector correlator for space- like 
q 2 was found for small temperatures (no OPE was used), 
a difference sum rule analysis yields an increase of the p - 
meson mass and a decrease of the mass of the a\ -meson 
with temperature. This was interpreted as the process of 
chiral symmetry restoration. The analysis of Ref. |l5| is 
meaningful for small temperatures (T < f % ) and for a T 
independent Lr consistent with the recent lattice simula- 
tion of Ref. (l9| indicating no sign of a T dependence up 
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to T « T c of the chiral transition. Our result for the T- 
evolution of the p -meson mass is consistent with the one 
obtained in the analysis of Ref . jD| . In addition, qualita- 
tive arguments based on the instanton model |30|] imply 
a cancellation of the effects of a decreasing constituent 
quark mass and a lowering of the instanton mediated at- 
traction between constituent quarks for the T-evolution of 
the p -meson mass. 

Depending on so(0), our result for the gluon conden- 
sate exhibits a 25-30% decrease at T w 160 MeV which 
is compatible with the lattice data for two quark flavor 
QCD of Ref. ||. The analysis showed, that the results are 
sensitive to the anomalous scaling of mass-dimension six 
operators under a change of the renormalization point (see 
Fig. 3). As Fig. 3 indicates the T-evolution of the gluon 
condensate is universal in a sense, that it is practically 
independent of the initial value so(0). 

There are still open questions concerning the inclusion 
of nonscalar, mixed operators at mass-dimension six in 
the thermal OPE. So far we must be content with the 
hope, that in the future their pionic matrix elements can 
be estimated from deep inelastic lepton scattering off the 
pion target Q). 

To conclude, our analysis indicates, that in the Gibbs 
average the T-dependence of the vacuum matrix elements 
of scalar operators, as introduced via a scaling relation in 
so, produces a sizable decrease of the gluon condensate 
and a moderate increase of the p -meson mass up to tem- 
peratures of about 160 MeV, where the analysis breaks 
down ( and the confidence in a good accuracy of the ap- 
proximations made is not high anymore). 
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"Struktur und Wechselwirkung von Hadronen und Kernen" (DFG 
Mu705/3), and T.G and A.F. were supported by the BMBF 
(contract no. 06 Tii 887). One of us (R.H.) would like to ex- 
press his gratitude towards Markus Eidemueller for reading the 
manuscript and suggesting some improvements. 



A 

Under the premises that T be considerably less than the 
fundamental QCD scale A, and that the renormalization 
scale Qo is taken of the order of 1 GeV we may argue in 
favor for the assumption of an implicit dependence of the 
scalar condensates on T via sq(T) in the chiral limit. 

At T — the most fundamental description of the 
QCD vacuum structure by the appearance of nonvanish- 
ing vacuum averages of local, scalar, and gauge invariant 
operators may in the chiral limit only involve the scale A 
and a renormalization scale Qo- With the number of inde- 
pendent parameters kept fixed, we hypothesize that there 
is a linear relation between A and y^so(0)y. In this pic- 
ture the constant of proportionality expresses this chan- 
nel dependence. In the following we will argue that for 

7 There is then no other, independent scale to allow for a dif- 
ferent behavior since so cannot depend on the external variable 
Qo. 



T > and provided that T is considerably smaller than A 
it is approximately possible to define a T-dependent fun- 
damental scale At such that the thermal vacuum averages 
of scalar operators have the same functional dependence 
on At and Q as in the case of T = 0, where they depend 
on A and Qo- Then there should be the identical linear re- 
lation between At and \J sq{T) as between A and \J s (0) 
at T = 0. 

For the thermal vacuum average of O renormalized at Qo 
we may write 

(0|0|0)r(Qo) = f(T/A, T/Q , A/Q ) x Q d , (50) 

where / denotes a dimensionless function of its dimension- 
less arguments, and d indicates the mass dimension of O. 
We define At by demanding 

/(0, 0, At/Qo) = f(T/A, T/Q , A/Q ) . (51) 

Solving Eq. (|l|) for A T yields 

A T = A T {A, Qo, T) = A T (T/A, T/Q , A/Q ) x A , (52) 

where At is a dimensionless function of its dimensionless 
arguments. Since in our applications the renormalization 
scale Qo is of the order of 1 GeV and with T considerably 
less than A, we can neglect the dependence of At on T/Qo 

A t (A, Qo, T) » A T (T/A, A/Q q ) x A , (53) 

where At again is a dimensionless function of its dimen- 
sionless arguments. Let us now expand At in powers of 
A/Qo 

A T = A x {c (T/A) + Cl (T/A)A/Q + ■■•)■ (54) 

Since At=o = A it follows that 

c (0) = 1 , a(0) = , (» > 1) . (55) 

Then the expansion of Eq. (|54|) should be well approxi- 
mated by its first term if T is considerably less than A, 
namely by 

A T ~ Ax c {T/A) . (56) 

Hence we obtain a well defined (since Qo-independent) 
T-dependent scale At which yields the same functional 
dependence of the thermal vacuum average on At and Qo 
as that of the T = vacuum average on A and Qq. 



B 

Here we derive the decomposition of Eq. (|2l|). For the 
operator C| 2 one simply uses 

qi a t a q = qLl a t a qL + q~Rlat a qR 

to obtain 

0| i2 = na s (P 4 + 2 P 2 ) . (57) 
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Using 

qjaj5t a q = q~Llat a qL ~ qRj a t a qR , 
yields the following decomposition for 0g 1 

O s 6>1 = na s [(^ Lla t a T 3 i; L ) 2 - 

2i; Lla t a T 3 ij L $ Rla t a T 3 ij R + 



(58) 



We are only interested in the flavor singlet part of the 
right-hand side, for which we also write 0g x 



O 



6,1 



1 

— na 
3 



[{i>Li a t a T b ii L ) 2 - 

2^Ll a t a T b ^ L ij R -f a t a T b ij R + 
^Rl a t a T b i) R ) 2 ] 
1 



-Tra s 
3 



(j Rla t a T b ij R ) 2 -2P 5 ] . 

For the sake of brevity we only consider the following op- 
erator 

Ol(r) = (MR)l a t a T b i; L{R) ) 2 . (60) 
Writing out color and flavor indices (in this order), using 

T ij T mn = 2(<5i„(5j m — -J^5ij5 mn ) , 

with N c — 3 and Nt — 2 one obtains 



O 



L(R) - 

4 [$L(R)Av7a8 lm S^\4> L ( R )jv$L(R)AKja8jl8 UK -lp L ( R ^ m x 
!^L(R),i^la8i m 8^L{R),]u'^L{R),lKla5jl8^\'ll}L{R),mX 
^L{R),i^a^ij^^\i l L{R)jv^L(R),lKlothm5vK'^L{R),mX 
^L(R),i^la5i ] 5^i) L i yR ^ 3v ^ L ^ R )^ K ^ a 8i m 8^x^ L ( R ^ m x 



Applying the Fierz transformation 1 14 



^L(R),lla^L(R),2^L{R)3la^L(R)A = 
^L(fl),l7aV'L(i?),4V'L(fl),37aV'L(fl),2 (62) 

to the first and the third line of Eq. (|rJl| ) and relabelling 
m <-> j in the second line, the sum of the second and the 
third line of Eq. ( |6l] ) reads 

5 - 

- Q'>pL(R),ila^L(R),j^L(R),jla^L(R),i , 

where the summation over flavor indices is implicit. This 
can easily be rewritten as 
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Putting everything together we finally obtain 
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6,1 



32 
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Figure 1: Temperature evolution of the p -meson mass. 
The solid lines correspond to the case (b) for So(Q) = 
1.2 GcV 2 (m p (T = 0) = 694 MeV), s (0) = 1.5 GeV 2 
(m p (T = 0) = 745 MeV), and s (0) = 1.8 GcV 2 (m p (T = 
0) = 776 MeV). For s (0) = 1.5 GeV 2 the case (o) is as- 
sociated with a dashed line. 



Figure 2: Temperature evolution of the pQCD thresh- 
old s - The solid lines correspond to the case (b). For 
so(0) = 1.5 GcV 2 the case (a) is associated with the 
dashed line. 



Figure 3: Temperature evolution of the gluon condensate 
normalized to its zero temperature value. For the case 
(b) we take s (0) = 1.2 GeV 2 , s (0) = 1.5 GeV 2 , and 
So(0) = 1.8 GcV 2 corresponding to dot-dashed, solid, and 
long dashed lines, respectively. For so(0) = 1.5 GeV 2 the 
case (a) is associated with a dotted line. 
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Fig. 3. 



